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Abstract. Motion of a point mass in gravitational fields of the Sun and of the galac- 
tic disk is studied. Fundamental features of the motion are found by investigating 
the time-averaged differential equations for orbital evolution. Several types of possi- 
ble orbits are mathematically exactly derived in a strictly analytical way. The relation 
Q,3 p2 _ j(^gp^ jg^ ^p) between semimajor axis a and period P of the change of osculating 
orbital elements is found (the index denotes initial values of the quantities). 

Due to conservation of energy in potential fields a is a constant. Moreover, the com- 
ponent of angular momentum perpendicular to the galactic plane is conserved. Due to 
these facts the system of equations reduces to two equations for either (e, lo), or (i, uj) 
(the length of the ascending node does not enter the equations for a, e, i, lo and is not 
solved here). 
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1. Introduction 



Motion of a point mass in gravitational fields of the Sun and of the galactic disk is 
important from the point of view of orbital evolution of cometary orbits. One of such 
numerical studies may be found, e. g., in Pretka and Dybczynski (1994). Since we want to 
find fundamental features of possible motions, we use time-averaged equations for orbital 



elements (e, i, lo) - conclusions may be applied to real situations for semimajor axes up 
to the order of lO'^ AU. 

2. Dynamical Model 

The dynamical model discussed in this paper is given by the following equations of mo- 
tion: 

y = — 3 y 

z = z — k z , (1) 

where iJ, = G Mq , k — A tt G g and g is the mean density of the mass of the galactic disk 
in the Solar neighbourhood. 

3. Time-averaged Equations for Orbital Elements 

The case /c = in Eqs. (1) corresponds to a Kcplcrian motion. If k is smaller than fi / r^ , 
we can consider the term /c z as a disturbing acceleration to two-body problem. The 
disturbing acceleration is given by 

a = (0, 0, - fc z) = (0, 0,- kr sin i sin 6) (2) 

in rectangular coordinates, where O = lu + /, / is the true anomaly. Radial, normal 
and transversal components of the disturbing acceleration are: 

an ~ a ■ eji — — k r (sini) (sinO) 

ajq — a ■ eiq ^ — k r (sinz) (cosi) (sinO) 

ttT — a ■ ex — — k r (sinz) (sinO) (cos0) . (3) 

Eqs. (3) enable us to obtain time-averaged evolutionary equations for orbital elements. 
The significant equations arc (fJ and T do not enter this system, and, in this paper, we 
are not interested in their evolution): 
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k ^ — (sinz) [sin (2 Cij)] e \/l 
V M 

< — > = k \ — [sin (2 i)] [sin (2 uj)] 



dt 8 y ^ vi - 



e^ 



duj 5 , /a^ I (sini)^ ^ e^ / • x2 1 A ;'l /,n 

<-r> = - ^ k J — { - — ^ (smw) - - \/l - e^} . (4) 



The first of Eqs. (4) accords with conservation of energy in the system. 
Eqs. (4) show that it is useful to define a dimensionless variable r; 




(5) 

Moreover, equations for eccentricity and inclination define the relation between these 
two quantities: 



Vl — e^ cosz= Jl — e^ cosio ? (6) 

where the index denotes initial quantities. Eq. (6) corresponds to conservation of 
the z— component - the component perpendicular to the galactic plane - of angular 
momentum together with the condition a = ag: dH / dt — rxa — rx ( —^r/r^ — k z e^r) 
= — r X e^r k z ^ H^ = Hn = constant. 

Eqs. (5) and (6) reduce Eqs. (4) to two equations: 

<-—> = - [sm[2u})\ (smz) 



dr 2 ^\ — e\ cosio 

1/5 — (sin a;) (l — Cq) (cosjq) + (sincj) (cosi) 



dhJ 

dr ^\ — eg (cosio) (cosi) 

4. Mathematical Treatment 



(7) 



Eqs. (7) form a complete set of differential equations - initial conditions ojo, io must 
be added, of course. We may interpret solutions of Eqs. (7) in terms of several types of 
orbits. 
At first 

n2 4 . 9\ / TT 3 7r\ di dio 

(cosio) = - (1 - e;^) A ( wo == - V cjo = — 1 =^ < ^ > = < -^ > == . (8) 

Eq. (8) define one type of orbits - stationary orbits. The remaining types correspond 
to periodic changes of the time-averaged elements. 

When considering periodic changes, one would await that tb may be positive as well 
as negative. One obtains, on the basis of the second of Eqs. (7) 

uj<0-^F=^- (1 - e^) < (sin uj sm if ll - ( -^) \ . (9) 

5 1^ Vsini/ J 

Eq. (9) states that if the case a; < may occur, the quantity lu can change periodically 
only in a small interval given by the values of e and i - oscillations in lu; the same holds 
for the quantity sin 6 = sin a; sini. This is the second type of orbits. (The consequence of 
Eq. (9) is e < sinz, and, (sin 6) > 1/5.) 



The final type of orbits corresponds to tlie case wlien tlie condition a; > is per- 
manently fulfilled: relation (9) can never hold during the orbital motion. In this case u 
increases monotonically from to 2 tt, & oscillates about zero. 

The last two cases should correspond to the types of orbits found in Pretka and 
Dybczynski (1994) in a numerical way - by numerical integration of Eqs. (1). The two 
types of possible orbits are separated by the conditions: i) always a; > holds ((sin 6) < 
1/5 <=> a; > 0, always); ii) oscillations in uj (oscillations in w 4^ (sin b) > 1/5). We sec that 
the results do not depend on semimajor axis and eccentricity (the statement in Pretka 
and Dybczynski (1994)), but on eccentricity e and e/sini. 

5. One consequence 

The consequence of Eq. (5) and periodical solutions of Eqs. (7) is that the relation a^ P^ = 
./(^oj^Oi^o) between semimajor axis a and period P of the change of the time-averaged 
orbital elements exists. This relation is interesting in comparison with the Kepler's third 
law: a'^ / r^ = /i / (4 TT^) ... T is the period of revolution in the ellipse. 

6. Conclusion 

We have found, on a firm mathematical basis, all types of orbits which exist in the outer 
part of the Solar System for a given model. Three types of orbits exist: i) always < lu >= 
(the same holds for other orbital elements), ii) < lu > is always positive, and, iii) < a; > 
oscillates during the orbital motion. We have completed and put into a correct form the 
"experimental" numerical results of Pretka and Dybczynski (1994), all in an analytical 
way. An interesting result a^ P^ = f{eQ,io,ujo) between semimajor axis a and period P 
of the change of time-averaged orbital elements was obtained. 
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